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Abstract 

We consider a molecular aggregate consisting of N identical monomers. Each monomer com- 
prises two electronic levels and a single harmonic mode. The monomers interact with each other 
via dipole-dipole forces. The monomer vibrational modes are bilinearly coupled to a bath of har- 
monic oscillators. This is a prototypical model for the description of coherent exciton transport, 
from quantum dots to photosynthetic antennae. We derive an exact quantum master equation 
for such systems. Computationally, the master equation may be useful for the testing of various 
approximations employed in theories of quantum transport. Physically, it offers a plausible expla- 
nation of the origins of long-lived coherent optical responses of molecular aggregates in dissipative 
environments. 
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I. INTRODUCTION 



Quantum systems coupled to a "bath" of harmonic oscillators are paradigmatic models 
in nonequilibrium quantum statistical mechanics llllel. The primary goal of the theory of 
quantum dissipative systems is the derivation of so-called master equations, which govern the 
dynamics of the reduced density matrix of a "system" of interest, starting from the Liouville 
- von Neumann equation for the combined system+bath ensemble. A well know example is 
Redfield theory (or modifications thereof f?!), which is valid in the limit of weak system-bath 



coupling [l 
coupling [9 



-|8[. A master equation can also be derived in the limit of the strong system-bath 



10|. More advanced treatments have lead to hierarchies of master equations 



which are valid beyond weak system-bath coupling and the Markovian approximation. 



With these methods, the time evolution of the reduced density matrix can be computed to 
high accuracy by the solution of truncated infinite chains of equations for auxiliary density 
matrices. 

There exist very few system-bath models for which "exact" closed-form quantum mas- 
ter equations have been derived. A master equation is called exact if the time evolution 
of all operators representing observables of the system is identical with the time evolution 
described by the Liouville - von Neumann equation for the system+bath ensemble. An ob- 
vious example is the so-called phase-noise case, in which the system Hamiltonian commutes 
with the system-bath coupling jlj]. Nontrivial exact quantum master equations are known 
for the free particle, the harmonic oscillator, as well as for two or identical oscillators 



bilinearly coupled to the harmonic bath 13Nl9|. The system-bath Hamiltonians that are 



quadratic in Bose or Fermi operators also allow for an analytical description of the system 
dynamics (see, e.g. 



20 



H22|). 



In the present communication, we present the derivation of an exact quantum master 
equation for a molecular aggregate consisting of A^ identical molecular monomers. Each 
monomer comprises two electronic levels and a single harmonic vibrational mode, which is 
coupled to a harmonic bath. This model is prototypical for the description of the transport 



of excitons, from quantum dots to photosynthetic antennae and DNA [4 .1231-12711. Apart from 



curiosity (as R. J. Baxter put it 



28 



"... the model is relevant and it can be solved, so why not 
to do so and see what it tells us?") our motivation is threefold. Conceptually, the derived 
master equation is, to our knowledge, the first exact master equation for a system with 
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electronic inter-monomer couplings, intra-monomer electron- vibrational couplings as well as 
vibrational dissipation. None of these couplings are assumed to be weak. Computationally, 
the exact master equation may be useful for the testing of various approximations (e.g., 
weak system-bath coupling, or weak electron-vibrational coupling, or weak inter-monomer 
coupling) which are frequently employed in theories of quantum transport 4. \24U2T\. |29|-|31|. 



Fundamentally, the master equation offers a plausible explanation of the origins of 



coherent optical responses of molecular aggregates in dissipative environments 
We use units in which ^ = 1. 
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ong-lived 

-kii- 



II. DERIVATION OF THE MASTER EQUATION 

Consider the total Hamiltonian H, which consists of the system (S) Hamiltonian, the 
bath (B) Hamiltonian, and their coupling, 

H = Hs + Hb + Hsb. (1) 

The system comprises N electronic two-level systems each of which possesses a vibrational 
mode: ^ 

Hs = H^ABl 5a', |X, - + E + + ^bIba) . (2) 

k=i ^ ^ 

Here, Bl and Ba' are the exciton creation and annihilation operators obeying the Pauli 
commutation rules [Ba,Bl,] = (5aa'(l — ^B^Ba')- Xk, Pk = —id/dXk, M, and Q denote the 
positions, momenta, masses, and frequencies of the harmonic oscillators. The parameter 
^ controls the strength of the intra-monomer exciton-vibrational coupling. The excitonic 
Hamiltonian ifg^. may depend parametrically on the distances \Xj — Xj/\. We need not 
specify Hex, but we require that it conserves the number of excitons and thus commutes 
with the number operator N: 

N 

[He,,N] = 0, N = Y,BlBk. (3) 

k=l 

We allow for variable couplings between the identical monomers, since the monomers are 
fixed (e.g., in protein matrices) and interact via (usually, dipole-dipole) forces, which depend 
on the relative positions and orientations of the monomers. 
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The system interacts with a harmonic bath via bihnear system-bath couphng 



H 



SB 



1=1 



Pi 

2mi 



+ 



N 
k=l 



niiUl 



(4) 



Here x,, Pi = —id/dxi, rrii, and Ui denote the positions, momenta, masses, and frequencies 
of the bath oscillators, and Cj are the system-bath coupling coefficients. All monomers 
are assumed to be identical, and the influence of the bath (jlj) on the system dynamics is 
determined by a single spectral density 



Nb 2 



2 ^ rriiUJ., 



CO 



(5) 



Otherwise, the masses M, frequencies Q, electronic couplings ^, and system bath-coupling 
coefficients acquire a fc-dependence, and an exact master equation cannot be derived. 
However, the master equation derived below remains true for different oscillators provided 
Mfe ~ Afc ~ ^fe ~ d^k for any k. 

In the description of energy-transport and relaxation phenomena in excitonic systems, 
the vibrational bath is usually assumed to couple electronic degrees of freedom directly, 
inducing fluctuations of the site energies 26|, |27|. In the present study, as well as, e.g., in 



Refs. y, 



42 



43| , each monomer is coupled to the bath through its vibrational reaction mode. 



As has been shown in 



44 



48], the two methods are, in principle, equivalent: By introducing 



an appropriate canonical transformation, one can switch from one description to another by 
incorporating the system modes into the bath or by singling-out several (high-frequency) 
modes from the bath and treating them explicitly. The present choice of the system-bath 
coupling through the monomer vibrational modes is motivated by three major reasons. 
First, it allows us to derive an exact master equation. Second, the explicit treatment of high 
frequency modes strongly coupled to electronic two-level systems allows us to assume that 
the remaining bath modes are coupled to the monomers rather weakly. This may become 
important for the extension of the present theory beyond the system-bath Hamiltonian @. 
Third, the explicit consideration of high-frequency vibrational modes facilitates the study 
of yibrationally- coherent effects in the system responses and spectroscopic signals (see, e.g.. 



42 



Let us introduce new variables specifying the system oscillators, the center-of-mass coor- 
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dinate R and the internal coordinates Qj 

1 ^ 

i?=-^X,-; Q,=X,+i-X,, j = l,2,...,iV-l. (6) 
The transformation ([6]) is described by the N x N matrix S: 

N N 

Qj = SjkXk, R = SNkXk- (7) 

k=l k=l 

The transformation back to the original coordinates Xk is given by the inverse matrix: 

N-l 

Xj = ''^{S~^)jkQk + iS^^)jNR- (8) 

k=l 

We were unable to derive an explicit analytical expression for for arbitrary N (for any 
finite N this can be done numerically). For the following derivations, it is sufficient to realize 
that S^^ obeys the conditions 

AT 

J2(S-'),k^0, {k = 1,2,..., N-l); {S-')jN^l, {j = l,2,...,N), (9) 

which are elementary consequences of Eq. ([8]). According to Eq. (jG]), the original momenta 
are connected to the new momenta Pr = —id/dR and Pqj = —id/dQj as follows: 

Pj = ^Pr-Pq, + Pq,-i- (10) 
The change of the system vibrational variables (|6]) is borrowed from the theory of Gaussian 



polymers (see, e.g., 49|). It is linear, non-singular (det S* = 1), but not canonical. The 
transformation ([6]) is just a technical tool, and we can return to the original canonical Xj 
representation at the end of the derivation. The choice (jH]) of the internal coordinates is 
natural for a linear array of monomers. Other choices may be preferable in different situa- 
tions. The transformation ([6]) can also be applied to the derivation of the master equation 
for N identical oscillators bilinearly coupled to a harmonic bath, generalizing thereby the 
approach of Ref . [lot . 

By insertion of Eqs. ([8]) and ffTOl) into Eqs. ([2]) and (jl]), we obtain the system Hamiltonians 
in the new variables 

Hs = Hi''^+HP, (11) 
(n) ( Pi MNQ'^R^ \ 



k=l k,l=l k=l 1=1 

(13) 

(by definition, Pq^q = Pq^n = 0). Here is the exciton number operator (j3]) and 

N 

Gki = Y.(S-')AS-%. (14) 
i=i 

The bath and the system-bath couphng Hamiltonians also split into contributions associated 
with the variables R and Q 

HsB + Hb = + Hf^i^ (15) 



where 



and 



Here 



i=l 



<) . G«^^. (17) 



Ab = 5^^, (18) 



.=1 ^^^^ 

SO that H^^j^ is the bath-induced renormalization of the potential energy (the so-called 
Lamb-shift). Combining the Hamiltonians, we have 

H = H^^^ +Hf^- (19) 

H^^^^Hf^+Hfl Hf^^Hf^ + Hf^ (20) 

(the subscript E indicates total system-|-bath Hamiltonians). Due to the requirement ([3]) 
the Hamiltonians (|T9|) commute 

[4^),4«)] = o. (21) 

Equations f|T9|) - fl2T|) summarize the first main result of the present paper. In words, they 

show that the parent Hamiltonian ([T]) can be transformed into the Hamiltonian fll9p . which 

is the sum of two mutually commuting Hamiltonians, h'^'^ and depends on the 

r(Q) 



internal coordinates Qj and is independent of the bath variables. The bath enters only 
through the bath- induced potential, . H^^^ depends on the center-of-mass coordinate 
R, through which the system is bilinearly coupled to the bath. The discrete degrees of 



freedom enter exclusively via the number operator A^. After the expansion in the 

eigenfunctions of A^, the Hamiltonian H-^ becomes harmonic. It is this latter fact which 
allows us to exactly integrate the bath out and derive an exact master equation. 

For doing that, let us introduce the total (system+bath) density matrix pit), which obeys 
the Liouville-von Neumann equation 

9,p(t) = -^[4^)+4'5\p(t)]. (22) 



18|, introduce a 



Since the Hamiltonians and commute, we can, following Ref. 
new density matrix 

9,p(«)(t) = -^[4^\p(^)(t)], (23) 
which is connected to the original density matrix via the unitary transformation 

p{t) = exp{-iH^h}p^^\t) exp{iH^^k}. (24) 

Since (as well as H^^) commutes with the number operator N, it can be expanded in 
the eigenfunctions of N: 

N 

H^,^^^Y.\^){n\H^l (25) 



n=0 

Here 

AR) - /^l rr(R) 



Ht/.u = (^1 4 1^) . N\n) =n \n) . (26) 
Therefore, Eq. ( l23l) can be rewritten in the form 

dtP^r^lit) = -tH^%p^il(t) + ^p^J^l{t)H^% (27) 



N 



(^)(t)^ \n){m\p^r^l{t). 



n,m=0 

Suppose that the density matrix at t = commutes with the number operator A^. Hence it 
can be expanded in the eigenfunctions of N as follows: 

N 

p(^)(0) = 5^|n)(n|p(f„)(0). (28) 

n=0 

Pnn{0) are not limited to pure states: they can represent a linear combination of pure states 
belonging to the same n. The initial condition (128|) is adequate for many practical purposes. 
Usually, the system is in its ground electronic state (n = 0) or can be promoted to its first 



excited electronic state (n = 1) by a laser pulse. Preparation of higher-order states (n > 1) 
is also possible via multiple and/or strong laser pulses. The level of description based on Eq. 
is appropriate for describing sequential (strong-pulse) spectroscopic signals, given the 



pulses are temporally well-separated 50|. If Eq. (^E^ is not fulfilled, we have to consider Eq. 
f l27p for n m. In this case, pmn{t) is propagated via the Hamiltonian H^l^ in the bra and 
via the Hamiltonian h'^^ ^ in the ket, and an exact master equation is impossible to 
derive. 

Given the initial condition ( 128|) . we restrict ourselves to the consideration of Eq. ( 127|) 
with n = m, 

dt~pi':.\t) = -z[Hgl~pi^\t)]. (29) 

The Hamiltonian H^l is the sum of the system Hamiltonian H'^!^^ f[T^ (in which me 
must substitute N by its eigenavalue n) and the system-bath Hamiltonian H^J^ ( JT6l) . Let 
us now define the renormalized system oscillator mass M = NM, the renormalized bath 
oscillator frequencies uf = Nuf, and the renormalized system-bath coupling coefficients 
Ci = Nci, as well as introduce the shifted center-of-mass coordinate R = R + nC,/{MQ'^), and 
the shifted bath coordinates Xi = Xi + n^c^ / {MVt^miUj]^) . After all these transformations, 
the Hamiltonian Hj^ ^ assumes the form 

Apparently, H^. is n-independent; it describes a harmonic oscillator of the mass M and 
frequency VL bilinearly coupled to the harmonic bath. The Liouville - von Neumann equation 
thus becomes 

5,p(^)(t) = -.[ijf ,p(^)(t)], (31) 



(t) being the total density matrix in the transformed variables. According to Refs. 



13 



16| , the bath can be integrated out in Eq. (13T]) exactly, yielding the master equation for the 



reduced density matrix p^^^t) = Ti B{p^^\t)}: 

dtp(^\t) = -t[Hf\p(^\t)] +M{t)p(^\t). (32) 

Here H^^^ = P^/(2M) -|- Mii^R^ /2 and the dissipation operator is defined as follows: 

^ -^a{t)[R\ p(^\t)] - ^b{t)[R, {Pn, p(^\t)}] 
+c{t)[R, [PR,p^'^\m - dit)[R, [R,p^'^\m. (33) 

8 



In the above equation, {., .} denotes the anti-commutator, and a(t), b(t), c(t), d(t) are real 
functions which are exphcitly defined in Isj^lGl; they are determined through the spectral 
density of the renormalized bath introduced in Eq. (15U]) : 

Nb -2 

«H-fE;55-*(^.-)-^«(;i) (34) 

i=l * 

{g{oj) is given by Eq. (j5])). Returning to the original variables R yields 

dtpi^it) = -^[Hitlpi^{t)] + 5?n(t)p(«)(t). (35) 

Here , Hg^^, and 3ft„(t) are obtained from p^^\ Hs^\ ^^'^ ^{t) by the replacement of R 
with /^ + </(M^]2). 

Combining Eqs. fl23|l and flMl) and using Eq. fl33|) . we obtain the desired exact master 
equation for the reduced (system) density matrix p{t) = Ti B{p{t)}: 

d,p{t) = -z[Hf^ + Hf\p{t)] + mpit), (36) 

3?(t) = Ylin 1^) ^nit) Equivalently, the master equation can be rewritten as 

dtpit) = -i[Hs + Hflpit)] + 3ft(t)p(t), (37) 

h'^^j^ (Eq. f|T7|) ) being the bath- induced renormalization of the system potential. Eqs. f p6|) 
and (1371) can be transformed back to the original Xj-representation, if desired. Very similar 
master equations can be derived assuming that the operators B\ and Ba' obey either the 
Bose or the Fermi commutation relations. 

From the computational point of view, H^^p cannot be factorized into commuting sub- 



Hamiltonians for the Qk- For small (see 5l|] for the vibronic trimer) a suitable matrix 



representation of h'^^j^ can be introduced. Alternatively, one can invoke the generalized 



Fulton- Gouterman transformation 



52^, 



53|, which diagonalizes h'^^j^ in the electronic Hilbert 
space, reducing the problem to the numerical evaluation of eigenvalues and eigenfunctions 
of the transformed Hamiltonian. 



III. DISCUSSION 



The master equation (137|) is mathematically exact and equivalent to the Liouville - von 
Neumann equation with the initial Hamiltonian ([T])-(jl]). According to the master equation 

9 



fl57]) . the density matrix in the Q subspace, p^^\t) = Tr(/j){p(t)}, does not experience 
dissipation. Indeed, 

V^)(t) = -^[4'3\p(«)(t)] (38) 

because the contributions due to Hg^'' and 3ft(t) are traced out (compare with Ref. 19|). 
The decoherence-free subspace [54 1 hence spans the entire Q subspace. Physically, this is a 
striking manifestation of quantum interference. Analogous effects have recently been studied 
in [55 1 . 

The presence of the decoherence-free Q subspace cannot be ubiquitous in real aggregates. 
The key simplifying assumption is that all monomers are coupled to a single harmonic bath 
(the Hamiltonian (j4])). More generally, we can assume that each monomer is coupled to its 
own bath, 

H,, + Hs = f:t(^ + '^ (H. -^J')- (39) 

^f— ^ \2mi 2 TUiuf J 

1=1 k=l 

The influence of these baths on the system dynamics is determined by the spectral density 
gjk{i^) ^ j,k < N enumerate different monomers). Apparently, the single bath (jl]) corre- 
sponds to Qjki^jj) = g{^) given by Eq. (|5]) and is thus equivalent to N fully correlated baths. 
Uncorrelated baths yield gjki^) = g{y>)5jk- Partially correlated baths can also be defined by 
introducing the correlation parameter < 7 < 1 as follows: gjk{^) = g{^){Sjk + 7(1 — ^jk))- 
The fully correlated and uncorrelated baths correspond to 7 = 1 and 0, respectively. If 
7 7^ 1, no decoherence-free subspaces exist and the entire system relaxes to equilibrium. 
The above analysis suggests that the relaxation in the R subspace has a characteristic rate 
u ~ 5'(w), while the relaxation in the Q subspace has a rate i/^ ~ z/(l — 7). 

It is impossible to derive an exact master equation in the case of partially correlated 
baths. However, the master equation (|37l) may be augmented with phenomenological dis- 
sipative operators which describe relaxation in the Q subspace. This can be done, for 
example, by switching to the eigenvalue representation in the Q subspace and introduc- 
ing the corresponding Redfield operator. Another option is to use the dissipation operator 
— z/^(l — pif ^Tr(Q){...}) 56|. Here Peg is the equilibrium Boltzmann distribution correspond- 
ing to the Hamiltonian 



The problem of bath correlations can be analyzed from a microscopic perspective [571] • If 
the harmonic potentials in the Hamiltonians Hsb and Hb are considered as linearizations 
of anharmonic interaction potentials between the particles of the system and the bath, one 
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arrives at the Hamiltonian because the identical particles of the system should interact 
via the same (in our case, harmonic) potentials with the bath particles. The domain of 
validity of such a primitive linearization of the interaction potentials is, however, limited 
to short times. On the other hand, we can consider the harmonic bath and the bilinear 



system-bath coupling in the spirit of a normal mode analysis [58|, |59|]. Each particle of the 
system then experiences different local potentials from the bath particles, no matter whether 
the particles are identical or not. This is tantamount to introducing different (possibly 
correlated) local harmonic baths for each monomer, which ensure relaxation of all {Q and 
R) degrees of freedom. On physical grounds, one expects that intramolecular vibrational 
baths should be close to the fully correlated limit, while environmental vibrational baths 
should be less correlated. 

Uncorrelated baths are the default choice in many simulations, although there exist strong 



experimental indications that this assumption is not universally applicable 32|-l34j. 



60 



6l|, 



and partially correlated baths are more appropriate for describing coherent energy and exci- 



ton transport |62h76| . The present analysis suggests that relaxation of molecular aggregates 
in the R and Q subspaces is governed by two different rates v and v^. For fully correlated 
baths, 1/^ = 0. One hence expects that < u for partially correlated baths. The slow relax- 
ation in the Q subspace manifests itself in slow decays of the populations and coherences of 
the aggregate density matrix. This may be one of the reasons of the existence of long-lived 



coherent optical responses in dissipative excitonic systems 321-141 1 
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